Conditions on the existence of localized excitations in 
nonlinear discrete systems. 
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1 Introduction 



Localization phenomena are of interest in nearly any branch of physics. In 
this paper we will deal with translationally invariant systems. In that spirit 
we deal with localization phenomena which appear due to intrinsic proper- 
ties of the underlying system in contrary to extrinsic sources (e.g. defects). 
There has been considerable success in the demonstration of the occurence 
of localized vibrations in Hamiltonian lattices The pure 

fact of the possibility of vibrational localization is astonishing because of the 
translational invariancy of the underlying lattice (i.e. no defects of any kind 
are necessary). The necessary localization condition has to be the nonlinear- 
ity of the lattice, since in the case of a linear lattice the problem is integrable 
and only extended degrees of freedom can be found. 

Most of the knowledge about vibrational localization is restricted to sim- 
ple one-dimensional lattices, usually with one degree of freedom per unit 
cell and nearest neighbour interaction. These systems belong to the class 
of Fermi- Pasta-Ulam (FPU) models (see e.g. 0). If one adds an external 
potential (field) which is periodic with the periodicity of the FPU system, 
one enters the world of Klein-Gordon lattices. The important property of 
vibrational localization in all these systems is that it is not of topological 
origin, i.e. we can consider a system with only one minimum in the potential 
energy function (the groundstate) and will find vibrational localization. If 
the potential energy function has several minima one can construct static 
kink solutions, i.e. static configurations of the system which link different 
groundstates (minima) of the model. Those static kink solutions can be ei- 
ther minima or saddle points in the potential energy function. The concept of 
vibrational localization is then extendable to those minima (kinks). One sim- 
ply considers excitations of the lattice above this (kink) minimum. Because 
the static kink himself breaks the translational invariance of the underlying 
lattice, a linearization and diagonalization of the fluctuations around kinks 
usually yields several eigenmodes which correspond to localized vibrations 
[§] , , ||10| , fill • The rest of the eigenmodes are deformed phonon-hke ex- 



tended degrees of freedom. Here we find a difference with the nontopological 
vibrational localization - namely localization in the linear limit. Of course it 
becomes only possible because the static kink already broke the translational 
symmetry. It is worthwhile to notice that the topological-induced vibrational 
localization is comparable to the well-known vibrational localization in a lin- 
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ear lattice with defects iT2 



Recently we were able to achieve progress in the understanding of non- 
topological vibrational localization |T^,[|n|. First we noticed that it is pos- 
sible to find nonlinear localized excitations (NLEs) which are essentially lo- 
cated on a very few particles. By analyzing the dynamics of the lattice on 
finite time scales we found that in general the NLE is a many-frequency 
excitation, where the number of frequencies n is equal to the number of par- 
ticles which are essentially involved in the excitation. By studying a reduced 
system which is given by the dynamics of the part of the lattice where the 
NLE is located (with proper defined boundary conditions) we showed that 
the NLEs in the original lattice correspond to regular motions of the reduced 
system on n-dimensional tori in its phase space. Moreover the NLE solution 
of the original lattice evolves on nearly the same torus. Thus it becomes 
possible to systematically study the NLE properties by checking the phase 
space structure of the reduced system. As a result we found that chaotic 
motion in the reduced system does not yield NLEs, as well as certain regular 
islands in its phase space which are well separated (by separatrices) from the 
NLE-regular islands. We were also able to attack the problem of movability 



of NLEs by considering a certain separatrix in the phase space |15]. Still 
the question remains, whether NLEs can be exact solutions of the lattice 
equations of motion. 

Much less is known about more complicated models. First we notice that 
the above cited approach in the simplest one- dimensional systems is easily 
extended to higher dimensions or to more than one degrees of freedom per 
unit cell. Indeed recent reports confirm the existence of NLEs in diatomic 
chains as well as in two-dimensional FPU systems 0,0. A 

systematic analysis of the NLE properties in two-dimensional Klein-Gordon 



lattices is carried out in |T^. Thus there is no doubt that by studying one- 
dimensional non-topological NLEs we are not restricting ourselves to exotic 
cases. The more complex reality will be covered. This is fundamentally 
different from the topological kink solutions. The reason is that the NLE 
solutions are not topologically induced. The only source of their existence is 
the nonlinearity of the underlying lattice. 

All presently known theoretical approaches to describe NLEs make the 
assumption that the NLEs exist. Then one can proceed in the description of 
their properties. In this paper we want to present an approach to the problem 
of the existence of NLEs. We will study the simplest one-dimensional cases. 
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As we have shown, this restriction is of no fundamental significance. We will 
prove several conditions when NLEs can not exist. Thus the remaining cases 
are the ones one has to choose if (possibly) NLEs exist. We will use the 
knowledge about the interpetation of NLEs in terms of actions and angles 
and consider a general ansatz. Then we reduce the problem of solving coupled 
ordinary differential equations to a (still highly complicated) set of coupled 
algebraic equations. The variables are certain Fourier-components. We will 
consider these algebraic equations as a mapping of a function defined on a n- 
dimensional lattice. By the definition of an NLE we show that the algebraic 
equations decouple in the tails of the supposed existing NLE solution. Then 
we analyse the decoupled mapping in the tails and calculate the eigenvalues 
of the linear map. We observe cases when periodic NLEs (n = 1) can not 
exist. We prove that strictly speaking there exist no solutions for n > 2, 
which implies that many-frequency NLEs are unstable. Still this bare fact 
allows for no strict conclusion about the typical decay time of many-frequency 
NLEs. We relate our studies to previous work on stability analysis of NLEs. 

2 Formulation of the problem 

We will consider the dynamics of a simple one- dimensional Hamiltonian lat- 
tice with the following Hamiltonian: 



Here Pi and Xi are canonically conjugated momentum and displacement of 
the /-th particle, V{z) and ^{z) are the potentials of the external field and 
nearest neighbour interaction respectively. We do not consider incommensu- 
rabilities between these two potentials, thus we assume that there exists at 
least one groundstate of |l| (minimum of the potential energy) such that with- 
out loss of generality X/ = for this groundstate. We specify the potential 
terms in ^ in form of an expansion around this groundstate: 




(1) 



OO -| 

V{z) = E -v,z' 

k=2 - 



(2) 



OO -| 



(3) 



k=2 
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The NLE solution is in its general form assumed to be given by the mo- 
tion of the phase space trajectory of on a n-dimensional torus |jl3[,[|l^ . 
Consequently the solution has to have the form: 



Xi{t) = E /«fcifc2...;i„e*^^^-'+'=^-^+-+'="-")* . (4) 

ki,k2,...,k„=—oo 

The localization property of ^ is defined by the boundary condition 

flkik2...k„ |/-*±oo^ . (5) 

Here we are excluding from the definition of an NLE localized pulses on 
carrier waves where the carrier wave doesn't decay far away from the center 
of the pulse. Since H is by assumption a solution of the equations of motion 

Xi = Pi = ~^ , (6) 

we can insert ^ into ^ and try to solve for the Fourier coefficients on the 
right-hand side of ^ Using 

Mt) = - E yd^)fik.k2..±y^''^'^''^'^---^'"^-^^ , (7) 

ki,k2,...,k„=—oo 

yj:{Q) = {kiUJi + k2UJ2 + ■■■ + kn^^nf , (8) 

k = {ki,k2,...,K) , U) = {uJi,UJ2, --^^n) , (9) 

we get Fourier series on the left and right hand sides of ^. The only possiblity 
of satisfying the obtained equation is to collect terms with equal exponents 
on both sides and to put the prefactors equal to each other. Then we obtain 
a highly complicated coupled set of algebraic equations for the Fourier coeffi- 
cients and the frequencies. Because we consider nearest neighbour interaction 
1^ we can formally write down the resulting set of equations: 

Mi^i,k = n{M,g},{M,_,g,}) . (10) 

Here we introduced a function ^ which is defined on a discrete n-dimensional 
lattice. The lattice is given by all combinations of {/ci, A;2; ^n} where each 
integer k^i varies from — oo to +oo. We have 

^Lk = flkik2...k„ ■ (11) 
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The mapping formally derived in [TO reminds us of the well-known two 



dimensional mappings which were used in order to study static kink proper- 
ties, commensurate-incommensurate transitions and breaking of analyticity 
pO[| . We can view ^ as a first step of implementing the fruitful ideas for static 
topologically induced structures in the dynamical problems of nontopological 
NLEs in Hamiltonian lattices. 

Let us study |10| in the tails of the NLE, i.e. for / ±00 where || holds 



by assumption. In the generic case f 2 and 02 from Q,^ will be nonzero. Then 
we can write down the mapping |l^ explicitely. We will do it without loss of 
generality for I — > +00. The corresponding formula for large negative / can 
be obtained by substituting /' = — /. We find 

M,-,irk = (^d^) + - ^i-iS ■ (12) 

Here we have introduced another function on the n-dimensional discrete space 
which is given by 

i^ki^) = 7 • (13) 

Equation [12| is linear and thus every component of M in the ra- dimensional 
discrete space decouples in this equation from all other components. Intro- 
ducing 

= M,r, - M,_,^^, (14) 

we finally arrive at a two-dimensional mapping for every component of ^ 
which reads 

Mi+iS = '^MMiS + GiS + ^iS ' (15) 

Gi^,:, = K^{mis + G,s ■ (16) 

Let us stress that under the assumption of an existing NLE solution the 
linearization of the map in the tails of the NLE is arbitrarily correct, if the 
distance from the NLE center is large enough. This mapping has a fixed 
point for ^ = and ^ = 0. It is characterized by the matrix A: 

We have 

det(A) = 1 . (18) 
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Thus the mapping p!5| , p!6| is symplectic and volume preserving. For the eigen- 
values of A we find 



A. = l + + , (19) 

A+A_ = 1 . (20) 



We can consider three cases: 

a) > : < A_ < 1 (21) 

i.e. A_ is real. Especially A_(Kg(u;) 0) — > 1 and A_(Kg(cJ) oo) — > 0. 

b) Kf:{Lu) < -4 : -1 < A+ < (22) 

i.e. A+ is real. Especially A+(Kg(cJ) — oo) and X+{kj^{uj) —4) —>■ 
-1. 

c) -4<K^(cu)<0 : |A+| = |A_| = 1 (23) 

i.e. A± are complex conjugated numbers on the unit circle. Consequently in 
cases a) and b) the fixed point of the mapping is a saddle point, i.e. there 
exists exactly one direction (eigenvector) in which the fixed point can be 
asymptotically reached after an infinite number of steps. In case c) the fixed 
point is a marginally stable elliptic point, i.e. starting from any direction the 
fixed point can be never reached after an infinite number of steps, instead 
the mapping will produce a (deformed) circle around the fixed point. Thus 
we find that case c) (|23|) contradicts the localization condition |^. 



3 Single-frequency localized excitations 

Let us consider n = 1. Then the NLE solution is periodic (cf. Equation 
p!3| can be simplified to 

The frequencies ujq for small-amplitude phonons around the considered ground- 
state of (where q is the wave number) are related to the parameters f 2 and 
<^2 by 

V2<^l<V2^ 402 . (25) 



Then it follows that case c) given in ^ is identical with 



c) kfujf = uj; . (26) 

We find that a single frequency NLE (periodic localized solution) can not ex- 
ist if any multiple of its fundamental frequency equals any phonon frequency. 
The reason is that we can not satisfy ^ and |^ simultaneously because of ^ 



In [|T3| ,|Ti| we have shown that under the assumption of the existence of 
a single frequency NLE its stability with respect to small amplitude phonon 
perturbations will depend on the fundamental NLE frequency. We found 
that if uJq/uJi = n/2, n = 0,1,2,... then the small perturbation will grow. 
Consequently the NLE would collapse. Here now we find that if n = 2ki 
(even n) then the NLE itself could not exist. U n = 2ki + 1 (odd n) then the 
NLE could exist, but would be unstable against phonon perturbations. 

One can interprete |2^ as a definition of nonexistence bands on the uji- 
frequency axis for different ki. Introducing a normalized frequency uji = 
uji/y^ and normalized interaction 02 = 4>2/v2 (only if V2 7^ 0) those bands 
are given by 



1 9 1 + 



For any finite 02 and small enough lji the nonexistence bands ^ will start 
to overlap. Consequently there will be always a lower bound on allowed 
NLE frequencies. Still some existence windows are possible in the phonon 
gap {uji < 1). However with increasing 02 more nonexistence bands will 
overlap and at the critical value 02 = 3/4 all nonexistence bands |2^ overlap. 
For all values of 02 > and any value of u)i > 1 + 402 condition ^ is 
not satisfied. That means that independent of the values of the parameters 
f 2 > and 02 > periodic NLE solutions are allowed with frequencies above 
the phonon band. From the above results it follows that model |l| always 
allows for periodic NLEs with frequencies above the phonon band. But if the 
model has a nonvanishing lower phonon band edge (^2 > 0) then periodic 
NLEs with frequencies in the phonon gap are allowed if the phonon band 
width is small enough compared to the lower phonon band edge. 

Now let us make some statements about symmetries of periodic NLEs 
if they exist. If the frequency of the NLE is above the phonon band then 
it follows < —4 for all ki. This corresponds to case b) in Then 

we have -1 < A+ < 0. Consequently -1 < M^^-^^/M^^ < for all / 
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in the tail of the NLE. Thus we find a coherent out-of-phase type of the 
motion of neighbouring particles in the tails of the NLE solution because of 
|: -1 < Xi+iit) / Xiit) < if defined. If the frequency of the periodic NLE 
is in the phonon gap {v2 > 0) things become more complicated. Namely 
there will always exist a certain finite integer kc such that for ki < kc it 



follows k,j:{lu) > which corresponds to case a) in |2T|. The corresponding 



Fourier components |TT| would yield in-phase type of motion in the tails of 
the NLE solution. However for all ki > kc the case b) in ^ applies. Those 
Fourier components would yield out-of-phase motion. Numerical findings 
indicate that usually the Fourier components decay very fast with increasing 
ki 1^. Then we could expect overall in-phase type of motion. However if 
the frequency Ui becomes smaller then it is well-known that the decay in the 
Fourier components with increasing ki slows down. Thus we have to expect 
a complicated mixture of in- and out-of-phase type of motion. 

4 Many- frequency localized excitations 



Let us consider n = 2. Then in analogy to 26 case c) in |2l applies if 



V2 < (kiLUi + k2UJ2f <V2 + 402 • (28) 

Now it is possible to show that there exists an infinite number of pairs of 
the integers (/ci,fc2) such that ^ is satisfied if the ratio uji/uj2 is irrational 
and f 2 > and f 2 > (cf. Appendix). Thus strictly speaking there exist 
no exact two-frequency NLEs. The proof for n > 3 is then straightforward 
and yields the same result. Of course this fact does not tell anything about 
decay times of many-frequency NLEs. It only states that many-frequency 
NLEs can not exist for infinite times. It seems to be logical to assume that 
the decay times are sensitive to the pair of the lowest integers {ki,k2) for 



which ^ holds in the case n = 2. Indeed numerical simulations ||13||,|T4| 
show extremely weak decay of two-frequency NLEs in Klein-Gordon chains, 
i.e. the characteristic decay time is several orders of magnitude larger than 
internal oscillation times. 
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5 Time-space separability 



Recently there were reports in the hterature where for systems of type |T] peri- 
odic NLE solutions with a property of time-space separability were proposed 
to exist ||21|. In more detail this property implies the existence of a master 
function G{t) in time such that the NLE solution can be given by 

Xi{t)=AiG{t) , Ai |/^±oo-^0 . (29) 

Without loss of generality one can set max(y4;) = 1. In terms of the intro- 
duced Fourier components in ^ ansatz imposes a rather strong symmetry 
on the Fourier components - namely they have to be equal to each other at 



different lattice sites up to a universal scaling number. If we insert ^ into 
the equations of motion || we find 

+ 00 

G{t) = -j:^kG'-\t) . (30) 

fc=2 

With 0,0 we can specify the constants rjk. 

Vk = ^[vkAf-' + MiAi - - {Ai+, - AO'-^}] . (31) 

If rik in ^ would depend on / then the differential equation for the master 
function G(t) in |3y would yield different solutions for different lattice sites 
(cf. e.g. 1^). Then we contradict the original ansatz Thus we have 



to assume that 1]^ in 31 is independent of /. Hence we generate a set of 



two-dimensional mappings for the vector {Ai, Ai_i) in pil for different k. 



Let us show that generically we can not satisfy all those mappings. For 
that we consider a Klein-Gordon lattice, i.e. 0^ = for k > 3, at least one 



Vk is nonzero for k > 3. Then we consider |31| for that specific k: 

rjk = VkA^-^ , k>3 . (32) 



From ^ it would follow that Ai = 1 for all / in contradiction to Thus we 
would have to conclude that either no NLEs exist in Klein-Gordon lattices 
or that they do not obey separability property The existence of periodic 
NLEs in various Klein-Gordon systems is verified numerically with very high 
probability ( i.e. currently periodic NLEs can be generated without any 
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measurable energy loss over 10^ time periods of the internal NLE oscillation, 
cf. Section 6). Consequently the ansatz ^is restricted to a certain subclass 
of systems which excludes the Klein-Gordon systems. Now let us consider any 
system such that for only (or at least) one integer k we have = (pk = 0. Let 
us assume that we have periodic NLE solutions which satisfy Perturbing 
this system with any small but finite ^ we can again consider the 
mapping for that k and yield Consequently no NLE solution would be 
allowed. Thus the separability ansatz p9| becomes a real exotic property. 



There is one case when we can at least expect that ^ holds. Namely 
when we have to satisfy only one mapping of the type This can hap- 
pen e.g. if Vk = for all k and 0^0 7^ for /eq < 4 and k^ even and 
0fc = for all k ^ ko. One can prove that if ^ is rewritten for that spe- 
cific case in terms of a mapping of the two-component vector (A;, Ai — A;„i) 
then the mapping will be nonconservative (it does not preserve phase vol- 
ume). A numerical test for ko = 4 yields that a periodic NLE solution can 
be found for rjkg = 9.5843773776314.... The amplitudes for that solution 
read ...,0.004795,-0.1657879,1,-1,0.1657879,-0.004795,.... First we no- 
tice that indeed for such special systems it can be possible to generate a 
NLE solution with time-space separability property |2^. However our result 
also indicates that even for such a highly nonlinear model (no linear disper- 
sion) the NLE if it exists has no compact structure, i.e. the amplitudes Ai 
are not exactly zero outside a finite volume of the solution. That fact can 
also be observed by looking at If outside of a finite volume of a solution 
all amplitudes would be equal to zero, by inverting the mapping we would 
generate zeros for all amplitudes in the excluded finite volume thus contra- 
dicting the ansatz that at least one amplitude is equal to 1. As a consequence 
the claimed compacton structure of a periodic NLE solution in ||2l[] is wrong 
and based on a simple calculation error. 

6 The decay of single-frequency NLEs in the 
tails 

Let us consider an existing periodic NLE for system |I] with f 2 > 0, 02 > 0. 
Then in the tail of the NLE (without restriction for large positive /) every 
Fourier component obeys a two-dimensional mapping given in |1^|1^. The 
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eigenvalue describing the decay is given by 0,0. Then the decay of every 
Fourier component will be governed by an exponential law ~ (sgn(A))'e''^'^''. 
The absolute value of A will depend on the number k of the Fourier compo- 
nent. For / — s> oo only the Fourier component with | A | closest to one will 
be present, i.e. all other Fourier components will decrease exponentially fast 
compared to the remaining one. The answer to the question which Fourier 
component 'survives' depends on the frequency uji and on the parameters V2 
and 02. Only if ui is above the phonon band then can we state that the 
'surviving' Fourier component is that with k = 1. 

In the following we will formulate two predictions for single-frequency 
NLEs in order to test them in a given realization. The first prediction was 
already formulated above - if an assumed periodic NLE (characterized by its 
fundamental frequency Ui) exists, the amplitude decay in the tails of the NLE 
(i.e. where the nonlinear terms in the equations of motion can be neglected) is 
governed by an exponential law with exponent /n|A|. Let us consider the case 
of a (j)^ chain, i.e. V{z) = — + 0.25z^ and ^{z) = O.bCz"^. Consequently 
we have V2 = 2, 02 = C. Periodic NLEs were investigated in |13[,[0 on finite 
time scales for the case C = 0.1. We use a particular realization with energy 
E = 0.256 (for details on the numerics the reader should consult the original 
work flSljliU)- '^^^ initial conditions on the lattice (3000 particles) are 
found from analyzing the elliptic fixed point of a reduced problem |T^. The 
frequency of the periodic localized object is found to be ui = 1.177 (cf. Fig. 10 
in [|14|). After a waiting time of T = 3 • 10^ the energy stored in five particles 
around the excitation center is still E{T) = 0.256. There is absolutely no 
drift (radiation) observable. The symmetric amplitude distribution in this 
assumed NLE solution is shown in Fig.l in a semilogarithmic plot (open 
squares). An evaluation of the eigenvalues corresponding to the different 
Fourier components (cf. p!9| - p^ ) reveals in this case that the eigenvalue with 
absolute value closest to one is given for k = 1 and reads A = 0.12465 
{ki = 6.1467). Consequently in the tail of the solution (where the nonlinear 
terms in the equations of motion can be neglected) we find an exponential 
amplitude decay with exponent logiQ|A| = —0.904. The dashed line is this 
prediction of the exponential decay. Note that in the semilogarithmic plot in 
Fig.l this decay law appears as a straight line. This line fits the measured 
amplitude decay down to amplitudes of 10"'^. The deviations for smaller 
amplitudes (farther away in the tail) are due to the fact, that the numerical 
realization of the NLE solution is always accompanied by the existence of 
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phonons. The phonons with nearly zero group velocities (wave numbers close 
to band edges) are practically not moving. These small phonon contributions 
increase (additively, as has to be expected in a linearizable tail of a NLE) the 
amplitudes of the particles in the NLE solution. We have no knowledge of any 
result in the literature, where from the measurement of the frequency of an 
assumed NLE solution the amplitude decay in the NLE tails was successfully 
predicted. 

The second prediction we wish to formulate is, that if the NLE frequency 
Ui is restricted to be in the nonzero gap of the phonon spectrum, then there 
exists a nonzero value of Ui^^ such that the exponential decay in the NLE 
tails will be weaker for all other frequencies Ui (still belonging to the gap). 
Let us explain why this statement follows out of the previous considerations. 
First the frequency Ui has to be larger than the phonon band width - else 
one (or more) of multiples of Ui will always lie in the phonon band. Secondly 
if Ui is slightly below the lower phonon band edge, then the decay will be 
very weak. Lowering the ui we increase the decay exponent, but since 2uji 
is coming closer to the upper phonon band edge, there will be a certain 
frequency when the decay in the tails will be governed by the second harmonic 
rather than the first harmonic of Ui. Let us calculate mm(l — |A(/tfc)|) with 
|A(Kfc)| < 1 and the parameters of the above introduced 0^ chain for different 
Ui- The result is shown in Fig. 2. We see indeed that there is a maximum 
in the decay exponent in this given example for uj^^^ = 0.938. As it can 
be observed from Fig. 10 in fl^, there are two NLE energies Ei = 0.336 
and E2 = 0.85 which correspond to this particular frequency uj["^\ The 
prediction is thus, that by varying the energy of the periodic NLEs, in a 
simulation one should observe maximum amplitude localization at these two 
energies Ei,E2. From the numerical runs reported in [|14| we can calculate a 
normalized amplitude entropy cxa which measures the amplitude distribution 
of the solution. It is defined in analogy to the energy entropy which is 
defined in |T^: 0"^ = — l/(ln(iV)) X^z ajna^ with ai = Ai/J2i^i- Here Ai 
is the amplitude of particle /. From the definition we have < cTq < 1. 
Maximum localization corresponds to a minimum of the normalized entropy. 
(Ja as a function of the NLE energy is shown in Fig. 3 (filled triangles and 
solid line). The predicted positions of the minima of (Ta{E) at Ei,E2 are 
given by vertical dashed lines. The predicted positions of the minima agree 
with the measured ones within the energy grid of the numerical experiment. 
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In order to demonstrate the usefulness of the above considered spatial 
decay laws, we have implemented a numerical method in order to solve the 
full nonlinear equations for the Fourier coefficients ^ for a periodic NLE 
(n = 1) on a finite lattice. We have choosen V2 = 1, = —1, ffc^2,4 = 0, 
02 = 0.1, 0fc^2 = and a system size of 100 particles with periodic boundary 
conditions. The NLE frequency was choosen to be below the phonon band: 
cui = 0.8. The resulting Fourier coefficients were computed on every lattice 
site with a maximum multiple of the fundamental frequency max = 30 
and a numerical error of lO"^''. The result is shown in Fig. 4 for thirty lattice 
sites around the NLE center in a semilogarithmic plot. Clearly a fc-dependent 
exponential decay of every Fourier component is found. Moreover in Fig. 5 we 
plot the measured exponents (slopes) as found from the numerical solution 
and compare them to the theoretical prediction. The agreement is very good. 

What happens if f 2 = and 02 = 0? We can give a particular answer for 
the special case 04 7^ and all other coefficients in being zero. This case 
was already discussed in the preceding section. If we know that a periodic 



NLE exists we can consider the mapping 31. We linearize the map near 



the saddle fixed point and find in leading order a decay of the amphtudes 
according to the law q-^-'^^p'^'- where a and b are positive numbers which 
depend on the parameters of the model Hamiltonian. 



7 Conclusion 

Nonlinear localized excitations (NLEs) might exist in nonlinear lattices for 
infinite times (i.e. they can be exact solutions of the equations of motion) if 
they are periodic in time and all multiples of the fundamental frequency are 
outside the phonon band. This is only possible because of the discreteness 
of the underlying lattice. A continuum system would have no upper phonon 
band edge, thus resonance would always be possible. Exceptions are classes 
of nongeneric systems where the resonance condition holds but the coupling 
between the NLE and the phonons is exactly zero. 

Let us note that we have not shown that periodic NLEs will exist on a 
given nonlinear lattice. We have restricted the possibilities of NLE solutions 
to time-periodic NLEs. 

If the NLE solution is described by two or more fundamental frequencies 
then there will be always resonance with phonons and thus those solutions 
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can not exist for infinite times, i.e. they are not exact solutions of the equa- 
tions of motion. Nothing can be said up to now about the hfetimes of such 
solutions. Numerical test show that lifetimes can be very large compared to 
internal oscillation times. If one considers nongeneric lattices without linear 
dispersion (phonons) then many-frequency NLEs might be exact solutions of 
the equations of motion. 

Space-time separability in the NLE solution can appear only in very non- 
generic cases. The construction of mappings for Fourier coefficients in higher 
dimensional Hamiltonian lattices seems to be more complicated but this is 
a technical question. The construction of the existence conditions for NLEs 
indicate that one can generalize them for higher dimensions too. 

Let us make some comments on the results rpesented in this paper. First 
we have assumed that the NLE solution (if it exists) is given by a regular 
motion on a torus in the phase space of the system (discrete spectrum). It 
is very hard to believe that NLE solutions can exist if the corresponding 
orbit belongs to a stochastic part of the system's phase space (continuous 
spectrum). Assuming the NLE solutions we search for are regular, all the 
subsequent steps in the presented analysis are free from any simplifications, 
approximations or conjectures. We have formulated the leading order am- 
plitude decay of an assumed NLE solution in its tails. We were thus able to 
formulate two nontrivial predictions and test them by comparing to numeri- 
cal experiments. 

Recently it was shown by using the properties of map [l^, that the ex- 
istence of NLE solutions is equivalent to finding a common point of two 



separatrix manifolds of the nonlinear map |2^. It was argued that the cor 



responding task is not overdetermined and should lead to a discrete set of 
solutions. Moreover in the case of an anharmonic Fermi-Pasta-Ulam (FPU) 
chain with homogenious potential Vk = 0, 02m = 1, <t'k^2m = (where m is 
any positive integer with m > 2) the existence of time-periodic NLEs could 
be strictly proved [P^. 

Since the presented analysis was carried out in the tails of the NLE so- 
lution (where the corresponding mappings can be linearized) the legitimate 
question arrises, why we do not have NLE solutions in a purely harmonic 
lattice. The reason for the nonexistence of a localized solution in such a 
linear lattice is that the two separatrix manifolds of the saddle fixed point 
never cross, because the eigenvectors of the fixed point uniquely define the 
positions of the two separatrix manifolds. In a nonlinear system, inverting 
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the discrete map (i.e. starting in the tail of an assumed existing NLE and 
iterating towards the center of the NLE) the increase of the amphtudes (or 
Fourier coefficients) will increase the contributions coming from the nonlin- 
earity. These nonlinearities will be the reason for the crossing of the separa- 
trix manifolds. We wish to emphasize that still these arguments can not be 
considered as a proof of the existence of NLEs. They provide the possible 
source of the NLE existence in the framework of the presented approach. 
The aim of the present work was to analyze the NLE properties in the tails 
assuming the NLE existence. The obtained results restrict possible NLE 
solutions to time-periodic solutions. 

Recently MacKay and Aubry have carried out a proof of existence of 
time-periodic NLEs for Hamiltonian networks of weakly coupled oscillators 



[ p^ . This strict proof remarkably incorporates Hamiltonian lattices with 
and without disorder. Consequently one can see, that (at least in the limit 
of weak coupling) the NLEs in a nonlinear lattice are continuously connected 
with localized modes of say a harmonic lattice with defects. These results 
should once and forever make it clear, that a harmonic lattice without defects 
is a very isolated system, and one should never use it as a reference system in 
order to judge results of nonlinear theories. In the present work the proof of 
MacKay and Aubry corresponds to the case Vk 7^ 0, 0^ <^ f2. According to 
[ p^ the spatial decay of NLEs is at least exponential, which can be specified 
to be strictly exponential in the present work without the requirement of weak 
coupling. Together with the strict proof of NLE existence in the mentioned 
anharmonic FPU chain |27] we thus have an analytical basis of the existence 
of NLEs in a broad class of systems. 

We were able to construct a mapping for the NLE solutions because of our 
knowledge about their proposed dynamical behaviour. This knowledge comes 
from the interpretation of the NLE solutions as motion of phase space points 
on tori in the phase space. The resulting mapping is an algebraic problem 
and can be considered as a task for itself. Here we meet numerous results 
for approximate NLE solutions which use the Rotating Wave Approximation 
(RWA) [|Tl,|2§],^,|§],[]g. Within the RWA all Fourier components higher 
than a certain hand-choosen order are set to zero. The resulting algebraic 
equations can be either solved selfconsistently or by means of useful iteration 
procedures. Only periodic NLEs were considered up to now. 

It is interesting to note that recently a combination of the lowest RWA 
order (taking into account only the lowest Fourier component) was refor- 
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mulated into a mapping [^. It would correspond to the mapping for the 
amphtudes of the NLE solution in the present work under the assumption of 
time-space separability. 

The fact that the existence of NLEs seem to be a generic property of (non- 
linear) Hamiltonian lattices implies that either experimental realizations of 
NLEs can be found with probability 1 or that the considered model classes 
are useless in describing reality. This promises some intriguing questions for 
the future. 
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A 

Let us prove that if a and h are two real numbers such that 

< a < 6 (33) 
and uJi and uj2 are two real numbers such that 

< Co-i < Co'2 (34) 

and the ratio uji/uj2 is irrational, then there always exists an infinite number 
of pairs of integers (/ci, ^2) such that 

a <| kiOJi + k2UJ2 \<b . (35) 

Let us introduce 

a = — , b = — , ^ = — . (36) 

UJ2 UJ2 UJ2 



Then ^ is aquivalent to 

a <\ kii + k2 |< h . (37) 
Let us choose the numbers c and d such that 

d < c < b , (i = min((c — a), (6 — c)) . (38) 
We can consider an arbitrary integer N such that 

N>—j, (39) 
nd 

where n = 1,2,3, ... and is fixed.. Then it follows that there exists at least 
one pair of integers (A;°, fcg) with 1 < A;^ < such that 

0</i<^ , ^=\k'^^ + k',\ (40) 

(cf. 1^). If we denote by m the integral part of a//i: 

m = [-] (41) 
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we find 



a < {m + n') ^ < b , 1 < n' < n . 



Thus our condition ^ can be always fullfilled if we clioose 
ki = {m + n')k° , k2 = {m + n')kl . 



(42) 



(43) 



Since tfiere is no restriction on tlie integer n in |3^ it follows from ^ that we 
generate an infinite number of pairs {ki, /C2) which satisfy our inequality 



19 



References 

[1] A. J. Sievers and S. Takeno. Phys. Rev. Lett. 61, 970 (1988). 

[2] J. B. Page. Phys. Rev. B41, 7835 (1990). 

[3] V. M. Burlakov, S. A. Kisilev, and V. N. Pyrkov. Phys. Rev. B 42, 4921 
(1990). 

[4] D. K. Campbell and M. Peyrard. in: CHAOS - Soviet American Per- 
spectives on Nonlinear Science, ed. by D. K. Campbell. American Insti- 
tute of Physics New York, 1990. 

[5] T. Dauxois, M. Peyrard, and C. R. Willis. Physica D 57, 267 (1992). 

[6] F. Fischer. Ann. Physik 2, 296 (1993). 

[7] J. Ford. Physics Reports 213, 271 (1992). 

[8] R. F. Dashen, B. Hasslacher, and A. Neveu. Phys. Rev. D 10, 4130 
(1974). 

[9] J. Goldstone and R. Jackiw. Phys. Rev. D 11, 1486 (1975). 

[10] Y. Wada and J. R. Schrieffer. Phys. Rev. B18, 3897 (1978). 

[11] R. Boesch and C. R. Willis. Phys. Rev. B39, 361 (1989). 

[12] H. Bottger. Principles of the Theory of Lattice Dynamics. Akademie- 
Verlag Berhn, 1983. 

[13] S. Flach and C. R. Willis. Phys. Lett. A 181, 232 (1993). 

[14] S. Flach, C. R. WiUis, and E. Olbrich. Phys. Rev. E 49, 836 (1994). 

[15] S. Flach and C. R. WiUis. Phys. Rev. Lett. 72, 1777 (1994). 

[16] M. Aoki. J. Phys. Soc. Japan 61, 3024 (1992). 

[17] S. R. Bickham, S. A. Kisilev, and A. J. Sievers. Phys. Rev. B47, 14206 
(1993). 



20 



[18] M. Aoki, S. Takeno, and A. J. Sievers. J. Phys. Soc. Japan 62, No. 12 
(1993). 

[19] S. Flach, K. Kladko, and C. R. Willis. Phys. Rev. E , in press (1994). 

[20] S. Aubry. Physica D 71, 196 (1994). 

[21] Yu. S. Kivshar. Phys. Rev. E 48, R43 (1993). 

[22] V. I. Arnol'd. Ordinary Differential Equations. MIT Press, Cambridge, 
1973. 

[23] S. Takeno and K. Hori. J. Phys. Soc. Japan 60, 947 (1991). 
[24] S. Takeno. J. Phys. Soc. Japan 61, 2821 (1992). 
[25] S. Wang. Phys. Lett. A 182, 105 (1993). 

[26] K. Chandrasekharan. Introduction to Analytic Number Theory. 
Springer- Verlag Berhn Heidelberg New York, 1968. 

[27] S. Flach. preprint (1994). 

[28] R. S. MacKay and S. Aubry. Warwick preprint 20/1994 (1994). 



21 



FIGURE CAPTIONS 



Figure 1 

Amplitude distribution for a localized state in the $^ chain (cf. text). 

Open squares - numerical result . 

Dashed hne - predicted decay law (cf. text). 



Figure 2 

Minimum of the distance of the absolute value of eigenvalue A from 1 as 
a function of the fundamental frequency Ui. Here |A| < 1, minimization is 
obtained with respect to all Fourier components for a chain with param- 
eters given in the text. 

Figure 3 

Normalized amplitude entropy Ua versus NLE energy for a chain (cf. 
text). 

Filled triangles - numerical result. 

Solid line connects the triangles and serves as a guide to the eye. 
Dashed vertical lines - predicted positions of minima of aa (cf. text). 

Figure 4 

Numerical solution for the Fourier components of the periodic NLE of a 
system with V2 = 1, V4 = —1, 02 = 0.1, N = 100, uji = 0.8 (cf. text). 
The absolute values of the components Aki are shown as functions of the 
lattice site I in a window of 30 lattice sites around the NLE center. The open 
squares are the actual results. The lines are guides to the eye and connect 
components with same Fourier order k. k increases from top to bottom as 
k = 1, 3, 5, 7, 23, 25. Fourier components for even k are zero because of 
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the symmetry of the potentiaL 
Figure 5 

Slopes of the hncs in Fig.4 as function of k (correspond to the exponents 
of the decay of the corresponding Fourier components, cf. text) are shown as 
open squares. The sohd hne connects the points of the theoretical prediction 
using the eigenvalues of the linearized map (cf . text) . 
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